Does hardcore interaction change absorbing type critical phenomena? 
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It has been generally believed that hardcore interaction is irrelevant to absorbing type critical 
phenomena because the particle density is so low near an absorbing phase transition. We study 
the effect of hardcore interaction on the N species branching annihilating random walks with two 
offspring and report that hardcore interaction drastically changes the absorbing type critical phe- 
nomena in a nontrivial way. Through Langevin equation type approach, we predict analytically the 
values of the scaling exponents, u± — 2, z — 2, a = 1/2, j3 — 2 in one dimension for all N > 1. Direct 
numerical simulations confirm our prediction. When the diffusion coefficients for different species 
are not identical, u± and (3 vary continuously with the ratios between the coefficients. 
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The study of nonequilibrium systems with trapped 
(absorbing) states has been very active in recent years 
jij. Models displaying absorbing phase transitions de- 
scribe a wide range of phenomena, in particular, epi- 
demic spreading, catalytic chemical reactions, and trans- 
port in disordered media. Furthermore, the absorbing 
transition is one of the simplest and natural extensions of 
the well-established equilibrium phase transition to non- 
equilibrium systems, which are still poorly understood. 

The concept of universality, which plays a key role in 
equilibrium critical phenomena, was shown to be appli- 
cable also to nonequilibrium absorbing transitions. Crit- 
ical behavior near an absorbing transition is determined 
by properties such as dimensionality and symmetry, and 
is not affected by details of the system. Finding a new 
universality class is difficult, and only a few classes of 
absorbing transitions are known jjj . 

Hardcore interaction between particles or kinks has 
been believed to be irrelevant to absorbing type critical 
phenomena, because the particle density is so low near 
an absorbing transition that the probability of multiple 
occupations at a site should be too small to be signifi- 
cant. This conventional belief leads to recent successes of 
field theoretical techniques using bosonic type operators 
[0-p|. However, it is well known that hardcore interaction 
does changes the asymptotic decay behavior of the parti- 
cle density in many multi-species diffusion-reaction mod- 
els near an annihilation fixed point Q. Since many ab- 
sorbing transition models can be mapped onto diffusion- 
reaction ones, it may seem natural to ask a question 
whether hardcore constraint changes the absorbing type 
universality classes in multi-species models. Despite re- 
cent efforts using fermionic formulation incorporating 
hardcore interactions ||J9], the effect of hardcore inter- 
actions is barely understood both analytically and nu- 
merically. 

In this Letter, we study the N species branching anni- 
hilating random walks with two offspring (7V-BAW(2)), 



introduced recently by Cardy and Tauber |q] . The model 
was solved exactly for all N > 1, using renormalization 
group techniques in bosonic type formulation which ig- 
nores hardcore interactions. We employ Langevin equa- 
tion type approach incorporating hardcore interactions 
and obtain analytically the values of critical exponents 
associated with the absorbing transition. It turns out 
that hardcore interaction drastically changes the univer- 
sality class in a nontrivial way and the critical exponents 
vary continuously with the ratio of diffusion constants of 
different species. Direct numerical simulations confirm 
our predictions. 

The 7V-BAW(2) model is a classical stochastic system 
consisting of N species of particles, Aj (i — 1, . . . , N). 
Each particle diffuses on a d-dimensional lattice with 
two competing dynamic processes: pair annihilation and 
branching. Pair annihilation is allowed only between 
identical particles (Aj + A; — > 0). In the branching pro- 
cess, a particle A$ creates two identical particles in its 
neighborhood (A, — > A; + 2 Aj ) , with rate a for i = j and 
rate a'/(N - 1) for i ^ j. 

For N = 1, this model exhibits an absorbing transition 
of directed Ising type (Zi symmetry) at finite branch- 
ing rate |l(] [l|]. The N species generalization imposes 
the permutation symmetry Sn between species. Like in 
the Potts type generalization of the absorbing transition 
models [0 , this model for N > 1 is always active except 
at the annihilation fixed point of zero branching rate. 

Critical properties near the annihilation fixed point 
have been explored exactly by Cardy and Tauber for 
N > 1 in the framework of bosonic field theory || . The 
upper critical dimension d c is 2. Using a perturbation 
expansion, they showed that the branching process as- 
sociated with a is irrelevant. For a = 0, it was found 
that the models for all TV > 1 are active for a' ^ and 
their scaling behavior near the annihilation fixed point 
{a' = a' c = ) forms a new universality class independent 
of N. For d < d c , the critical behavior is characterized 
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by the exponents 

is ± = l/d, z = 2, a = d/2, /3 = 1. 
Here, the exponents are defined as 
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where A = a' — a' c , £ is the correlation length, r the char- 
acteristic time, p(t) the particle density at time t, and p s 
the steady-state particle density. 

Even in the presence of hardcore interaction, the scal- 
ing exponents a and z should follow from the simple ran- 
dom walk exponents; z = 2 and a = d/z for d < d c ||. 
Near the annihilation fixed point, elementary scaling the- 
ory ensures (3 = v±za, which leads to /3 = v\d. We de- 
termine the value of v±_ through Langevin equation type 
approach. 

The particle density can change by branching processes 
and pair annihilation processes. If we start with a config- 
uration of very low particle density, the particle density 
initially grows by branching processes, (Ai — > Ai + 2Aj). 
In this growth regime, a newly created pair of offspring 
and its parent are far more likely to annihilate against 
each other than other particles in the system. Dynamics 
of such three particle configurations or "triplet" governs 
the growth behavior of the particle density and the inter- 
triplet interactions can be ignored. The particle density 
growth will be finally capped by pair annihilations pro- 
cesses of independent particles and the system reaches a 
steady state. 

We focus on the growth regime dominated by triplet 
dynamics, from which v± can be evaluated We only 
consider the case of a = 0, where a newly created pair is 
always dissimilar to its parent. The survival probability 
S(t) of the triplet of the same species decays much faster 
(~ t- 3 / 2 ) than that of different species, so the branch- 
ing process associated with a is irrelevant. Near a' = 0, 
the time evolution of the particle density of i-th species, 
Pi(t), is written as 



dpi 
dt 
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Pj(t) 



Lij(t - t')pj(t')dt' 



(3) 



where Lij(t — t')dt is the probability that an i-th species 
pair created by an j-th species particle at time t' , annihi- 
lates in an interval between t and t + dt. The two terms 
in the right-hand side represent the creation and annihi- 
lation process of a triplet, respectively. Pair annihilation 
contribution from independent particles is 0(p 2 ), which 
is ignored in the growth regime. 

The kernel Lij(t) is simply related to the survival 
probability Sij(t) of the triplet (Aj + 2AA as Lij(t) = 
—dSij(t) I 'dt. To keep the lowest order of a' in Eq. (^), 
we evaluate S%j(i) at a 1 = 0. When hardcore interaction 



is not present, a pair of Ai's does not see its parent Aj, 
so annihilate each other freely by random walks. In that 
case, it is well known that S(t) = S, 
totically as S(t) ^ 
becomes finite (S 



ij(t) decays asymp- 
t~ 5 with 5 = 1 - d/2 for d < 2 and 
= 0) for d > 2, irrespective of their 



diffusion constants |15|. However, with hardcore interac- 
tion, the pair annihilation process changes significantly 
due to an effective bias in their diffusive behavior, gen- 
erated by their parent particle Aj. The motion of Ai 
near Aj picks up a convective component with velocity 
proportional to t~ x l 2 , so the convective displacement is 
of the same order of diffusive displacement t 1 ! 2 . In this 
case, the competition between the convection and diffu- 
sion becomes nontrivial and the scaling exponent 5 de- 
pends continuously on the parameters of the system jl6) . 

We calculate the survival probability S(t) of a triplet in 
one dimension. With hardcore interaction, S(t) depends 
crucially on where to create two offspring with respect 
to their parent. When two offspring are divided by their 
parent (static branching) they have no chance to 

meet each other. The survival probability never decays 
(5 = 0). When two offspring are placed both to the left 
or both to the right side of the parent particle with equal 
probability (dynamic branching) fl^| , S(t) decays with a 
nontrivial scaling exponent. 

Consider three random walkers on a line, labeled as A, 
B and C . A is a parent particle that created two off- 
spring, B and C, to the right side of A. Two offspring B 
and C are of the same species, which is different from its 
parent A. Hardcore repulsion is present between A and 
B. B and C annihilate instantaneously upon collision. 
The calculation of S(t) belongs to the class of problems 
known as "capture process" fl1^|l7j| . 

Let the coordinates of the walkers be xa,xb and Xc, 
and their diffusion coefficients Da,Db and Dc, respec- 
tively. In our case, Db = Dc- It is useful to in- 
troduce the scaled coordinates = Xi/\fD~i, where 
i = A, B, C. Then we can map this triplet system to 
a single walker system with isotropic diffusion in three 
dimensional space (vAiVBiVc) Jl5|, 18|] . The walker sur- 
vives inside the wedge bounded by two planes: a reflect- 
ing plane P r of ' aVa — VDbVb and an absorbing 
plane P a of \JUbVb = VDcVc- 

The survival probability S(t) of an isotropic random 
walker in a d dimensional cone with absorbing bound- 
ary is known |19| ], In particular, S(t) in a wedge with an 
opening angle 9 asymptotically decays as t^/ 20 pg|| . In 
our case, one of the boundary planes, P r , is not absorb- 
ing but reflecting. The probability of the walker at P r is 
nonzero and there is no net flux across this plane. Us- 
ing this fact, one can easily show that our system should 
be equivalent to the system in a wedge bounded by two 
absorbing planes with twice bigger opening angle. 

We find that the survival probability of the triplet de- 
cays with the exponent 
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where O is an opening angle of the wedge and r = 
Da/Db- The exponent S monotonically decreases from 
1 to 1/2 as the diffusivity ratio r varies from to oo. At 
r = 1 (the same diffusivity for all walkers), S = 3/4. 

First, we consider the case that diffusion coefficients 
are identical for all species. The N coupled Langevin 
equations, Eq. (^|), can be simplified in terms of the total 
particle density, p(t) = ^Pi(i), as 

^ = 2<r'p{t) - 2a' [ L(t-t')p(t')dt', (5) 
dt Jo 

where L(t) = Lij(t) is independent of i and j. Taking 
Laplace transformation, we find 

sp(s) - p(0) = 2a'(l - L(s))p( S ) = 2a'sS(s)p(s), (6) 

where p(s) = J Q °° p(t)e~ st dt, and similarly L(s) and S(s) 
are the Laplace transform of L(t) and S(t), respectively. 
With S(t) ~ t- s , one can show that S(s) ~ s 5 ' 1 for 
6 > 0. 

The function p(s) has a pole in the positive real axis 
(T /i/(i-<5)_ ^hen the initial density p(0) is small, 
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the density p(t) increases exponentially as exp(s Q i). Us- 
ing the definition of the characteristic time r (Eq. (||)), 
we find 



1/S 



,-1/(1-5) 



(7) 



With (5 = 3/4 for the dynamic branching model, we ar- 
rive at v±z = 1/(1— S) = 4. Therefore we predict that the 
critical exponents for the dynamic branching A-BAW(2) 
model with hardcore interaction in one dimension are 



= 2,z = 2, a = l/2,/3 = 2, 



(8) 



which should be valid for all N > 1. For the static 
branching 7V-BAW(2) mo del, S = and v± = (3 = 1/2. 
Without hardcore interactions, branching methods do 
not matter and 8 = 1/2, which leads to Eq. (Q). 

We check the above predictions for the A r -BAW(2) 
model by direct numerical simulations for N — 2, 3 and 
4. We start with a pair of particles. With probability 
p, a randomly chosen particle (Ai) creates two offspring 
(2Aj) on two nearest neighboring sites (dynamic/static 
branching). The branching probability p is distributed 
as 7p for i = j and (1 — j)p/(N — 1) for i ^ j. Oth- 
erwise, the particle hops to a nearest neighboring site. 
Two particles of the same species at a site annihilate 
instantaneously. In case of models with hardcore inter- 
actions, branching/hopping attempts are rejected when 
two particles of different species try to occupy the same 
site. Critical probability p c = for all models considered 
here. 



We measure the total particle density p s in the steady 
state, averaged over 5 x 10 2 ~ 5 x 10 4 independent sam- 
ples for several values of A = p — p c (0.001 ~ 0.05) and 
lattice size L (2 5 - 2 11 ). We set 7 = 1/2. Using the 
finite-size scaling theory |2Q] 



(9) 



the value of v± is determined by "collapsing" data of p s 
with (3/v±_ = 1 (Fig. 1). Numerical data show that v± 
does not depend on N in all models as expected. Wc 
find = 1.00(5) for models without hardcore interac- 
tions, which agrees with the result by Cardy and Tauber 
With hardcore interactions, we find v± — 1.9(1) for 
the dynamic branching models and v± = 0.50(3) for the 
static branching models, which confirm our predictions 
within statistical errors. 

When the diffusion coefficients are not identical for dif- 
ferent species, Sij (t) decays with the exponent 6 depend- 
ing on diffusivity ratio r = Dj/D{. Instead of a single 
Langevin equation, we are then forced to deal with the 
N coupled Langevin equations. The solution of the sys- 
tem of equations is difficult in general, but the equations 
become quite simple for N = 2. 

Laplace-transformed coupled equations for N — 2 be- 
come 



spi(s) - pi{0) = 2a'sSi2 {s)p 2 {s), 
sp 2 (s) - p 2 (0) = 2a'sS 2 i(s)pi(s). 



(10) 



We take p 2 (0) = as an initial condition and solve the 
equations for pi(s): 



spi(s) - pi(0) = ^a' 2 sS 12 S 2 iPi{s). 



(11) 



Note that S\ 2 (t) decays with exponent S(r) with r = 
D 2 /D 1 and S 21 (t) with 5(l/r), see Eq. (|). From the 
pole position of pi(s), we arrive at 



ux(r) = 



1 



2-6(r) -S(l/r) 



(12) 



The exponent 6(r) ranges from 1/2 to 1, but S(r) + S(l/r) 
varies only slightly with r. It ranges from 3/2 to 1.5255, 
so v±(r) varies only within a few percent. Due to rather 
large statistical errors (~ 10%), we could not confirm nu- 
merically the r dependence of i>i_. However, it is clear 
from our derivation that v± should vary continuously 
with diffusivity ratio. Although we were not able to ob- 
tain a similar expression for vj_ for general N, we expect 
that v± varies continuously but only slightly with r for 
all N > 1. 

In summary, we showed that hardcore interaction 
in the 7V-BAW(2) model changes its universality class 
in a nontrivial way. Details of branching methods 
(static/dynamic branching) and also the diffusivity ratios 
between different species change drastically the absorb- 
ing type critical phenomena. We find that, for all N > 1, 
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the dynamic branching models with hardcore interaction 
form a new universality class, different from the mod- 
els without hardcore interaction. Especially, the scaling 
exponents vary continuously with the diffusivity ratios. 
The static branching models with hardcore interaction 
form yet another new universality class. Numerical sim- 
ulations confirm most of our predictions, but large scale 
simulations are necessary to measure the diffusivity ratio 
dependence of the scaling exponents. 

The present analytic method to study the effect of 
hardcore interaction can be applied to a wide range of 
multi-species diffusion-reaction models near the annihi- 
lation fixed point. Our analysis implies that many multi- 
species models with hardcore interaction may exhibit a 
nontrivial absorbing phase transition with continuously 
varying exponents. 
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FIG. 1. Data collapse of p s L Plv± - against AL 1 /" 1 with p/u± = 1 for various system size L = 2 5 , . . . , 2 11 for N = 2 BAW(2) 
models (a) without hardcore interaction, (b) with hardcore interaction (dynamic branching) , and (c) with hardcore interaction 
(static branching). Best collapses are achieved with (a) u± = 1.00(5), (b) 1.9(1), and (c) 0.50(3), respectively. 
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